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We report results of spin diffusion measurements in normal phase of liquid 3 He confined in nafen. Nafen is 
a new type of aerogel and it consists of AI2O3 strands which are nearly parallel to one another at macroscopic 
distances. We examine two samples of nafen with different porosities using spin echo techniques. Spin diffusion 
of 3 He along and across the strands was measured. The aerogel alignment is clearly evident from observed 
spin diffusion anisotropy. A theory describing this effect is developed and compared with the experiment. 


1. INTRODUCTION 


2. THEORY 


Superfluid 3 He in high porosity aerogel is a model 
system to investigate the influence of impurities on un¬ 
conventional superfluidity. Silica aerogels which consist 
of a nearly chaotic array of Si 02 strands are used in most 
of such experiments. An important parameter for theo¬ 
retical models of superfluid 3 He in aerogel is the mean 
free path (A) of Fermi-liquid quasiparticles which can be 
determined from measurements of spin diffusion coeffi¬ 
cient ( D ). At high temperatures (T > 20 mK) the den¬ 
sity of quasiparticles is large, so A and D follow the bulk 
Fernri-liquid behavior, i.e. A oc T~ 2 and D oc T~ 2 . At 
low enough temperatures (T < lOmK) aerogel strands 
limit the mean free path and the spin diffusion coeffi¬ 
cient, so that at T ~ 1 mK the density of quasiparticles 
is so small that values of A and D are fully determined 
by the array of aerogel strands and do not depend on T. 
Such behavior was observed in the first measurements 
of spin diffusion in 3 He confined in nearly isotropic 95% 
and 98% open silica aerogels [U [2]. In some recent ex¬ 
periments with superfluid 3 He in aerogel a new type 
of aerogel was used ®m®. The remarkable feature of 
this aerogel, called “nematically ordered” (N-aerogel), is 
that its strands are oriented along the same direction. 
There are two types of N-aerogel: “Obninsk aerogel” 
produced by Leypunsky Institute (Obninsk, Russia) [B, 
which consists of AlOOH strands and nafen (7) which 
consists of AI 2 O 3 strands. In the limit of T = 0 the 
strong global anisotropy of N-aerogel should result in 
anisotropy of 3 He spin diffusion. For example, in “Ob¬ 
ninsk aerogel” with overall density ~ 30 mg/cm 3 the 
spin diffusion along the strands is about twice as fast as 
that in perpendicular direction [5], 

Here we present results of the spin diffusion measure¬ 
ments in 3 He confined in nafen, which is much denser 
than “Obninsk aerogel” and compare the obtained re¬ 
sults with the theory we developed. 
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An approximate description for the weak field spin 
diffusion in bulk 3 He was first given by D. Hone [5] 


vf A ; 

T 


■1 VFA n^p°-\ dM 
Jm =- + 


(1) 


where M is the magnetization, jm is respective current, 
vf is the Fermi velocity, Fq is the Landau Fermi-liquid 
parameter, and A is a quasiparticle mean free path. Dif¬ 
fusion in anisotropic aerogel is described by a general¬ 
ized equation 
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where D lm is the spin diffusion tensor. It is important 
to note, that no mean free path can characterize diffu¬ 
sion in anisotropic medium. Diffusive flux m in such 
system depends on the gradient direction and they are 
not necessarily parallel to each other. On the other side, 
a mean free path is a scalar, it is by its nature averaged 
over quasiparticle distribution which is isotropic even in 
anisotropic aerogel. 

Suppose that the stationary kinetic equation 
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holds separately for each individual spin component. In 
the left hand side of the kinetic equation © a local equi¬ 
librium function should be substituted. This gives [9; 
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The magnetization here and below is taken in magneton 
units 7 h = 1 . 

Calculation of D lm is greatly simplified at low tem¬ 
perature so that the collisions between Fermi quasipar¬ 
ticles can be neglected. This means that the scattering 
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of quasiparticles on the aerogel alone is responsible for 
the collision integral 

/[n] = J dp' (w( p, p'K(l - n) - w( p', p)n(l - n ')), 

where n = n(p), n' = n(p'), and ic(p,p') = 7c(p',p) is 
the scattering p' —>• p probability. Since the scattering 
by aerogel strands conserves the spin and the energy, 
the collision integral vanishes for equilibrium distribu¬ 
tion function of the true quasiparticle energy rio(e). It is 
therefore possible to linearize the collision integral using 
Sh = n — no(e) 

I[n] = J dp'ur(p, p') (Sh 1 — Sh). 

Furthermore, the solution of the kinetic equation can be 
sought for in the form 

Sh=^ X (p), P = P /p, 

where x depends on the direction of p only. The kinetic 
equation becomes 

- = VF J dcr (P’ p'Xx' - X), (4) 

where der is the differential scattering cross section per 
unit volume. 

The flux jm also vanishes for no(e) and can be lin¬ 
earized as follows 

= / jSyw 6 " = -(Hy / <5) 

We now have all tools ready for the diffusion tensor 
calculation. The differential cross section der depends 
on the microscopic aerogel structure and the properties 
of quasiparticle scattering by aerogel strands. Let us 
represent the aerogel as an array of infinite cylindrical 
strands with diameter d oriented in z direction and ran¬ 
domly distributed in xy plane with the surface density 
N. The porosity of such structure is therefore 


7 rNd 2 



For the scattering on the cylinder walls we consider two 
opposite limits: diffuse and specular reflection. 

21 Diffuse Reflection 

After diffuse reflection all information about the ve¬ 
locity direction of the incident particle is lost. The scat¬ 
tering cross section on a unit wall area is 


where the unit vector m is the outer normal to the sur¬ 
face element. For a system of z-aligned cylinders, this 
vector is uniformly distributed in xy plane. The differ¬ 
ential scattering cross section for this system is obtained 
by integration 


da = ^ j dmiS(m z ) dcr m = 


Nd 
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sinf? sin<?' I sin A — A cos A| dp', 


where A = </> — <!>' £ (—7r, 7r), and the spherical angles 
(0,<f>) and (O',$') correspond to p and p' respectively. 

Generally, the diffusion tensor D lm has three princi¬ 
pal values. Due to the axial symmetry of the system, it 
can be characterized by the mere two distinct compo¬ 
nents ZXII and D^. It is therefore sufficient to calculate 
the diffusion in two directions: along and across the 
aerogel axis. In the former case the kinetic equation Q 
is 

- ip cos 0 = J(x - x) dtr. (7) 

If we substitute © in © and use the fact that the dis¬ 
tribution function does not depend on the polar angle 
</>, we get 
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The solution of this equation 

cos 0 
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should be substituted in © to get the flux 
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Similar procedure can be used to investigate the lat¬ 
eral diffusion. If ip m has only x component, then the 
kinetic equation has the form 


—ip sin 0 cos<f) 



X) der. 


dcr m 


To, (mp) > 0; 

< 0, (mp') < 0; 
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It has a solution 
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which leads to the magnetization flux 


3. DETAILS OF EXPERIMENTS 
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The diffusion anisotropy is described by the ratio 
D l l/D^ = (tt 2 + 16)/8 « 3.23 . 


(9) 
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22 Specular Reflection 

Elementary scattering cross section on a unit smooth 
surface is 


dcr m = - dp'(pm) 


0 , 


(mp) > 0; 


8 (p' — p + 2m(pm)) otherwise. 


For the array of cylinders we again integrate over pos¬ 
sible m orientations 


dcr = I dm S(m z ) da m = 


= — 6 ( 0 - 9 ') 
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The delta function guarantees the conservation of the 
^ component of momentum. This means that strictly 
parallel specular cylinders do not obstruct the spin 
flow along the system axis. Corresponding component 
of the diffusion tensor (infinite in our model) is 
limited either by strand irregularities or quasiparticle- 
quasiparticle collisions. Here we formally calculate only 
lateral (finite) component of the diffusion tensor. 

We again write down the kinetic equation 

fix' -X) to. 


It has a solution 


X = 


3 ip 
4Nd 


COS (j), 


We have used two samples of nafen with overall 
densities 90 and 243mg/cm 3 (nafen-90 and nafen-243). 
Their structure was investigated in m- The samples 
consist of strands with diameter of « 8 nm and « 9 nm 
correspondingly. The density of bulk AI 2 O 3 is 4g/cm 3 , 
this gives for the porosities of the samples 97.8% (nafen- 
90) and 93.9% (nafen-243). 

The experimental chamber used in the present work 
was similar to the chamber described in [3]. It is made 
of Stycast-1266 epoxy resin and has two separate cells. 
The cuboid samples with characteristic sizes of 4 mm 
are placed freely in the cells, so that ~ 70% of each cell 
is filled with nafen. 

Experiments were carried out using spin-echo tech¬ 
nique in magnetic field « 280 Oe (corresponding NMR 
frequency is « 900 kHz) at pressure of 2.9 bar. In or¬ 
der to avoid a paramagnetic signal from solid 3 He on 
the surface of strands, the samples were preplated by 
~ 2.5 atomic layers of 4 He. The strands of the samples 
were oriented parallel to the external steady magnetic 
field H. Two gradient coils were used to apply the field 
gradient in directions parallel and perpendicular to the 
strands. The necessary temperatures were obtained by 
nuclear demagnetization cryostat and were measured by 
a quartz tuning fork. The temperature was determined 
in assumption that the resonance linewidth of the fork 
in normal 3 He is inversely proportional to the temper¬ 
ature El. We calibrated the fork at high temperature 
where the diffusion coefficient in aerogel should be the 
same as in bulk 3 He (using [ 12 ] as a reference). 

Spin echo decay curves were obtained by standard 
two-pulse method: we measured amplitude of the echo 
after 7r/2 — t — 7t/2 pulses, where r is the delay between 
pulses. The measurements were carried out at temper¬ 
atures from 1.4 mK up to 60 mK for two directions of 
magnetic field gradient (parallel and perpendicular to 
the direction of aerogel strands) and at several values of 
the gradients (0.24 4- 1.25 0e/cm). 


4. RESULTS 

Expression for the spin echo amplitude 


which results in the magnetization flux 
3 t/> p% 


Jm ~ 


where 


4 Nd (27Hi) 3 


J d(j)d0 sin 2 6 cos 2 4> = 


d£ = (i + fs) 


37 tvf 37t 2 (1 + Fq ) vf d 


32 Nd 


128 


1 -p 


• (12) 


I = Iq exp(—2 t/T 2 — At 3 ) (13) 

can be found from Bloch-Torrey equations [TJ. With 
an obvious generalization for an anisotropic media, the 
coefficient A here is given by 

A = - 7 2 D lm G l G m , (14) 
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x = 2/3 (yG 3 )Y (s/cm 2 ) 

Fig. 1: Spin echo decay in 3 He in nafen-90 for dif¬ 
ferent field gradients applied in direction perpendicular 
to aerogel strands. G x = 0.70 0e/cm (open circles), 
0.53Oe/cm (filled circles), 0.37Oe/cm (open squares), 
0.24Oe/cm (open triangles). H = 278 Oe. T ~ 2.0 mK. 
Solid line is the best fit of the data at x > 100 by 
Ea. (1131) . Inset: Typical echo signal of 3 He in nafen-90. 

T ~ 2.0 mK, G x = 0.53 Oe/cm, r = 10 ms. 

where 7 is the gyromagnetic ratio, G l is the magnetic 
field gradient. 

Typical echo signal of 3 He in nafen is shown in the in¬ 
set in Fig HI In order to determine the value of spin dif¬ 
fusion coefficient, spin echo amplitudes should be mea¬ 
sured for different r and then fitted by Eq. (USD- In this 
procedure the term with X 2 can be neglected, because 
observed relationship between I/I 0 and G 2 r 3 does not 
depend on field gradient at all used temperatures. An 
example data set is shown in FigJT] Note that at x < 80 
experimental points in FigJT] deviate from the linear de¬ 
pendence. It is due to the presence of bulk 3 He out¬ 
side the aerogel sample. At low temperatures the spin 
diffusion in bulk 3 He is greater than that in aerogel. 
However, for the same reason the relative contribution 
of bulk 3 He into the total echo signal rapidly decreases 
with the increase of r, so we determined the value of 
spin diffusion coefficient D(T ) of 3 He in aerogel from 
the data at relatively large x where they follow the lin¬ 
ear dependence (x > 100 for FigO}. 

The measured temperature dependencies D(T ) for 
two orientations of the gradient are shown in Fig [2] for 
nafen-90 and in Fig|3]for nafen-243. In order to obtain 
a value of spin diffusion coefficient in zero temperature 
limit (D = D( 0)) these dependencies were fitted by the 
equation 

(15) 



T (mK) 

Fig. 2: Temperature dependence of the spin diffusion 
tensor: Dg 0 (T) (open circles) and Dg 0 (T) (filled trian¬ 
gles) in nafen-90. 



T (mK) 

Fig. 3: Temperature dependence of the spin diffusion 
tensor: D\ i3 {T) (open circles) and D^ i3 (T) (filled tri¬ 
angles) in nafen-243. 

where Df, oc T~ 2 is the diffusion coefficient in bulk 3 He, 
which is determined only by collisions between quasi¬ 
particles. Solid lines in these graphs are fits by Eq. (fl5l) . 
dashed lines - diffusion coefficient in bulk 3 He (extrap¬ 
olation to P = 2.9 bar using the data presented in l2l). 
Obtained principal values of the spin diffusion tensor 

are: 

• in nafen-90 

Ug 0 ss 0.049 cm 2 /s, 

Dg 0 ss 0.015 cm 2 /s, 

Dl/Di 0 « 3.3. 


D- 1 (T) = D^ 1 (T) + D~ 1 , 
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• in nafen-243 

-D 243 « 0.029 cm 2 /s, 

D 243 ~ 0.0036 cm 2 /s, 
^243/^243 ~ 8 - 1 - 

We estimate the accuracy of these values as ±10%. 


5. DISCUSSION 

In isotropic system the spin diffusion coefficient is 
(see Eq.([T])) 

D= l ^v F (l + F 0 a ). (16) 

As explained in Section 2 no mean free path can prop¬ 
erly characterize the case of globally anisotropic aero¬ 
gel. However, to account for the experimental data it is 
convenient to introduce zero-temperature effective mean 
free paths A|| and Ax using Ea- dHll) as the definition. If 
v F = 5397 cm/s and F 0 = —0.717 [H] then for nafen-90 
we get Ay ps 960 nm, Ax ~ 290 nm and for nafen-243 
Ay « 570 nm, Ax ~ 70 nm. 

In our experiments the strands of nafen were pre¬ 
plated with 4 He. In this case the reflection of quasi¬ 
particles on the surface is believed to be specular or at 
least partly specular (15 16 . For specular reflection 
is limited by imperfect strands alignment and irreg¬ 
ularities on their surface. Therefore we can not directly 
compare the measured ratio k = D H /D 3 - with proposed 
theory for specular case. However, for nafen-243 we ob¬ 
tain k « 8.1 which is much greater than it is expected 
for diffuse one (see Ea. (fIUl) l. This means that the re¬ 
flection in nafen-243 is predominantly specular. 

The values of Dj 3 and Dg can be calculated from 
Eqs. ©, dm using values of p and d for our samples. 
These theoretical values should weakly depend on vari¬ 
ations in orientations of the strands and can be com¬ 
pared with experiment. For the case of nafen-90 we 
obtain Dg ss 0.013 cm 2 /s and Dp ss 0.010 cm 2 /s while 
the experimental value Dg 0 ss 0.015 cm 2 /s. For nafen- 
243 Dg ss 0.0052 cm 2 /s, Dp ss 0.0043 cm 2 /s and the 
experimental value 3 ~ 0.0036 cm 2 /s. 

We can conclude that the experimental values quali¬ 
tatively agree with the theoretical model. The quantita¬ 
tive difference may be due to experimental errors in de¬ 
termination of D, errors in nafen parameters, and vari¬ 
ations of strands diameters (by ± 20 % as follows from 
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